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Abstract. Recent results from the Pierre Auger Observatory show energy dependent chemical 

\a ■ composition of ultrahigh-energy cosmic rays (UHECR) with a growing fraction of heavy elements 

at high energies. This points to a non-negligible contribution of the Galactic sources, such as past 

GRBs and other rare but powerful stellar explosions in the Milky Way. The effects of diffusion in 

the Galactic magnetic fields alter the observed composition and render the flux of UHECR isotropic, 

up to a few per cent anisotropy in the direction of the Galactic Center, as well as some small-scale 

anisotropy with "hot spots" due to the locations of the most recent/closest bursts. 
i-Q ■ 
Q_i' Keywords: cosmic rays, gamma ray bursts, galactic magnetic fields 

PACS: 98.70.Sa, 98.70.Rz, 98.35.Eg 

Composition of ultrahigh-energy cosmic rays (UHECR) measured by the Pierre 
Auger Observatory (PAO) shows a steady increase of the mean nuclear mass with en- 
,_i \ ergy between 2 EeV and 35 EeV [1, 2]. This unexpected result is difficult to reconcile 

>• | with the usually assumed extragalactic origin of UHECR. The lack of plausible sources 

in the Milky Way and the lack of Galactocentric anisotropy of the arrival directions of 
UHECR are usually presented as evidence for extragalactic origin of UHECR above 
10 18 eV. However, it was recently shown that, if the cosmic rays of 10 18 — 10 19 eV are 
nuclei produced in the Milky Way, the effects of diffusion in turbulent Galactic micro- 
Gauss magnetic fields can explain both the change in composition and the approximate 
anisotropy [3]. As for the plausible sources, there is a growing evidence that long GRBs 
are caused by a relatively rare type of supernovae, while the short GRBs probably re- 
sult from the coalescence of neutron stars and/or black holes. Compact star mergers 
undoubtedly take place in the Milky Way, and therefore short GRBs should occur in our 
Galaxy. Although there is some correlation of long GRBs with star-forming metal-poor 
galaxies [4], many long GRBs are observed in high-metallicity galaxies as well [5, 6, 7], 
and therefore one expects that long GRBs should occur in the Milky Way. Less powerful 
hypernovae, too weak to produce a GRB, but can still accelerate UHECR [8], with a sub- 
stantial fraction of nuclei [9, 10]. If the observed cosmic rays originate past explosions 
in our own Galaxy, PAO results have a straightforward explanation [3]. 

GRBs have been proposed as the sources of extragalactic UHECR [11, 12, 10], 
and they have also been considered as possible Galactic sources [13, 14, 15]. It is 
believed that GRBs, hypernovae, or other stellar events capable of producing UHECR 
could have happened in the Milky Way at the rate of one per toRB ~ 10 4 — 10 5 years 
[16, 17, 18, 19]. Such events have been linked to the observations of INTEGRAL, Fermi 
and PAMELA [20, 21, 22, 23]. As illustrated in Fig. 1, diffusion depends on rigidity, 




FIGURE 1. For each species, there is a critical energy Eq i for which the Larmor radius Rj is equal 
to the magnetic coherence length l c . For E <C Eqj, the mean free path of the diffusing particle is I ~ Iq, 
and D{(E) = l c /3. For E 3> Eqi, the particle is deflected only by a small angle 9 ~ h/Ri, and, after A: 
deflections, the mean deflection angle squared is 6 2 ~ Ic(Iq/Rj) 2 . The corresponding diffusion coefficient 

isD,(£)-(4) 2 , for£»£ ,,. 



and, therefore, the observed composition can be altered by diffusion [24, 3]. 

Changes in composition due to a magnetic fields have been discussed in connection 
with the spectral "knee" [24], and also for a transient source [25]. A simple model of 
diffusion [3] shows how the diffusion coefficient affects the observed spectrum of the 

species "?'". If all species are produced with the same spectrum n- = tiq <* E~Y at 
the source, their observed spectra are altered by the energy dependent diffusion and by 

the trapping in the Galactic fields, so that, instead of nf rc ' = hq, one obtains 

Hence, the composition becomes energy dependent. Indeed, at critical energy Eq^, which 
is different for each nucleus, the solution (1) changes from oc E~ r to oc E~ r ~ 2 because 
of the change in Dt(E), as discussed in the caption of Fig. 1. Since the change occurs at 
a rigidity-dependent critical energy £0.; = eEoZt, the larger nuclei lag behind the lighter 
nuclei in terms of the critical energy and the change in slope. If protons dominate for 
E < Eq, their flux drops dramatically for E > Eq, and the heavier nuclei dominate the 
flux. The higher Z ; , the higher is the energy at which the species experiences a drop in 
flux. 

One can also understand the change in composition by considering the time of diffu- 
sion across the halo is tj ~ R 2 /Dj. The longer the particle remains in the halo, the higher 
is the probability of its detection. At higher energies, the magnetic field's ability to delay 
the passage of the particle diminishes, and the density of such particles drops precipi- 
tously for E > Eqi. Since £ ; is proportional to the electric charge, the drop in the flux 
occurs at different energies for different species. 

The model of Ref. [3] gives a qualitative description of the data. To reproduce the 
data more accurately, it must be improved. First, one should use a more realistic source 
population model. Second, one should include the coherent component of the Galactic 
magnetic field. Third, one should not assume that UHECR comprise only two types of 
particles, and one should include a realistic distribution of nuclei. Finally, one should 
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FIGURE 2. UHECR spectra according to Calvez et al. (shown here for values of parameters that 
differ from those in Ref. [3]) from Galactic sources assumed to produce 75% protons and 25% iron, 



with identical spectra ^ E~ . The source distribution traces the distribution of stars in the Galaxy. The 
magnetic field was assumed to be ~ 10/iG, coherent over Iq — 100 pc domains. The overall power and the 
iron fraction were adjusted to fit Pierre Auger data points [26] (shown). 



include the extragalactic component of UHECR produced by distant sources, such as 
active galactic nuclei (AGN) and GRBs (outside the Milky Way). A recent realization 
that very high energy gamma rays observed by Cherenkov telescopes from distant 
blazars are likely to be secondary photons produced in cosmic ray interactions along 
the line of sight lends further support to the assumption that cosmic rays are copiously 
produced in AGN jets [27, 28]. For energies £>3x 10 19 eV, the energy losses due to 
photodisintegration, pion production, pair production and interactions with interstellar 
medium become important and must be included. The propagation distance in the 
Galaxy exceeds 10 Mpc, so that the Galactic component should exhibit an analog of 
GZK suppression in the spectrum. The extragalactic propagation can also affect the 
composition around 10 18 eV [29]. 

Galactocentric anisotropy for a source distribution that traces the stellar counts in the 
Milky Way is small [3] . Although the anisotropy in protons is large at high energies, their 
contribution to the total flux is small, so the total anisotropy was found to be < 10%, 
consistent with the observations. The latest GRBs do not introduce a large degree of 
anisotropy, as it would be in the case of UHE protons, but they can create "hot spots" 
and clusters of events. 

Our model [3] makes an interesting prediction for the highest-energy cosmic rays. 
Just as the protons of the highest energies escape from our Galaxy, they should escape 
from the host galaxies of remote sources, such as AGNs. Therefore, UHECR with E > 
3 x 10 19 eV should correlate with the extragalactic sources. Moreover, these UHECR 
should be protons, not heavy nuclei, since the nuclei are trapped in the host galaxies. 
If and when the data will allow one to determine composition on a case-by-case basis, 
one can separate £>3x 10 19 eV events into protons and nuclei and observe that the 
protons correlate with the nearby AGN. This prediction is one of the non-trivial tests of 



our model: at the highest energies the proton fraction should exist and should correlate 
with known astrophysical sources. 
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